Tensor-force-driven Jahn- Teller effect and shape transitions in exotic Si isotopes 



Yutaka Utsuno, 1 Takaharu Otsuka, 2 ' 3,4 B. Alex Brown, 4,5 
Michio Honma, 6 Takahiro Mizusaki, 7 and Noritaka Shimizu 3 

Advanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan 
2 Department of Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
3 Center for Nuclear Study, University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing MI 48824, USA 
J Department of Physics, Michigan State University, East Lansing MI 48824, USA 
6 Center for Mathematical Sciences, University of Aizu, 
Ikki-machi, Aizu-Wakamatsu, Fukushima 965-8580, Japan 
7 Institute for Natural Sciences, Senshu University, Tokyo, 101-8425, Japan 
(Dated: January 20, 2012) 

We show how the shape evolution of the neutron-rich exotic Si and S isotopes can be understood as 
a Jahn- Teller effect that comes in part from the tensor-driven evolution of single-particle energies. 
The detailed calculations we present are in excellent agreement with known experimental data, 
and we point out of new features that should be explored in new experiments. Potential energy 
surfaces are used to understand the shape evolutions. The sub-shell closed nucleus, 42 Si, is shown 
to be a perfect example of a strongly oblate shape instead of a sphere through a robust Jahn- Teller 
mechanism. The distribution of spectroscopic factors measured by 48 Ca(e,e'p) experiment is shown 
to be well described, providing a unique test on the tensor-driven shell evolution. 

PACS numbers: 21.60.Cs,21.30.Fe,27.40.+z,21.10.Pc 



Among the new frontiers of nuclear physics, one of the 
most important is the evolution of single-particle energies 
in nuclei far from stability that led to dramatic changes 
in the location of magic numbers and types of collectivity 
relative to nuclei near stability [lj . One of the ingredients 
is the tensor-force driven changes in single-particle ener- 
gies leading to tensor- force driven shell evolution @, Q , 
with many experimental examples, e.g., 

We shall show in this Letter that tensor-driven shell 
evolution plays a critical role in the rapid shape change 
as a function of neutron and/or proton number, including 
triaxial and 7-unstable shapes. In particular, we show 
for the first time [B[ how this shape change at low energy 
can be related to the Jahn- Teller effect p, 0] , where the 
degenerate intrinsic states led to a geometric distortion 
that removes the degeneracy and lowers the energy. 

We study the structure of Si and S isotopes as ex- 
amples, using shell- model calculation with SDPF-MU 
Hamiltonian to be described. The sd — pf cross-shell 
proton- neutron interaction is important, and is given by 
a monopole-based universal interaction (Vmv) P( except 
for minor details. We show a test of this Hamiltonian 
by the spectroscopy of K isotopes. We then demonstrate 
that we can describe also the distribution of spectroscopic 
factors measured in the (e,e'p) experiment j8|, exhibiting 
the first theoretical reproduction of them. We then show 
levels, B(E2) values and potential energy surfaces (PES) 
of Si and S isotopes for even A=22~28 in comparison 
with experiments [§-[l5| , following an intuitive picture of 
the robust appearance of the oblate deformation in 42 Si. 
We shall show that not only the property of 42 Si but also 
the rapid change of structure of all these Si and S isotopes 
are naturally described. 



We first outline the details of present shell-model cal- 
culations. The sd and pf shells are taken as the valence 
shell, whereas no excitation between them is included 
for the sake of simplicity. The interactions within each 
of these shells are based on existing interactions: USD 
0] (GXPF1B [13]) for the sd {pf) shell. Regarding 
the monopole interaction [H, , as pointed out in Ref. 
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must be corrected from its USD value. 
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We made a shift, AV T d =% d5/2 = -3 X AU C 
-0.7 MeV, consistently with V M u and SDPF-M [2C 
The monopole- and quadrupole-pairing matrix elements 
(O/7/2O/7/2JVI 0f 7/2 0f 7/2 )j=o,2 are replaced with those 
of KB3 [18( for a better description of isotopes of N ~ 
20. The cross-shell part, most essential for exotic nuclei 
of interest but rather undetermined so far, is given basi- 
cally by Umu of Q with small refinements stated below. 
The tensor-force component of Vmu is nothing but the 
7r + p meson exchange force. It has been used in many 
cases [1, [3] , and has been accounted for microscopically 
under the new concept of Renormalization Persistency 
[2ll | . The central-force component is fine tuned from Vmu 
with density dependence similar to the one in [22| , so that 
its monopole part becomes closer to that of GXPF1 [23j |. 
We include the two-body spin-orbit force of the M3Y 
interaction [3]. Following USD and GXPF1B, all of 
the two-body matrix elements are scaled by A~ 3 . The 
single-particle energies (SPE) of the sd shell are taken 
from USD, and those of the pf shell are determined by 
requesting their effective SPEs on top of 40 Ca closed shell 
equal to the single-particle energies of GXPF1B. 

The Hamiltonian thus fixed is referred to as SDPF-MU 
hereafter, and is diagonalized by the mshell64 code [HI]. 
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FIG. 1: (Color online) Evolution of 1/2^ level of K isotopes 
measured from 3/2+ level. Experimental data (filled circles) 
are compared with theoretical results including (solid line) or 
excluding (dashed line) the cross-shell tensor force. 



Note that none of the cross-shell monopole interactions 
are fitted directly to experiment, in contrast to other re- 
cent interactions 2(| 27 1. 

We first study proton single-particle states of odd- A 
K (Z=19) isotopes. Figure Q] exhibits 1/2 J" levels for 
7V=20~28. If valence neutrons occupy the lowest possi- 
ble orbits (filling configuration), the 1/2+ level relative 
to 3/2+ corresponds to the 0g? 3 / 2 -1si/2 gap. Since this 
is the case with N = 20 and is almost so with N = 28, 
the lowering of the 1/2+ level from N = 20 to 28 re- 
flects considerable reduction of this gap. Experimentally 
observed reduction (~3 MeV) is reproduced remarkably 
well by SDPF-MU Hamiltonian. This reduction is due 
to the proton-neutron monopole interaction Vqj 7 a od 3/2 



that is more attractive than VJf ( ™ , 

0/7/2, J-Sl/2 



the difference is 
0.44 MeV (at A = 42), out of which the tensor and the 
central forces contribute, respectively, by 0.21 MeV and 
0.22 MeV. The central force yields a stronger attraction 
between 0/7 / 2 and 0d 3 /o because of the similarity of their 
radial wave functions [3j . 

We move on to the splitting between the 0e? 3 / 2 and 
0c?5/2 single-particle energies. The 0ri5/2 single-particle 
strength is highly fragmented due to its large excitation 
energy (6-7 MeV). Spectroscopic factors of the proton sd- 
shell orbits have been measured for 48 Ca by one-proton 
removal through (e, e'p) reaction Q. The left panel of 
Fig. [2] displays the experimental values in comparison 
to those obtained by the present calculation, where the 
usual overall quenching factor 0.7 is used [28| . The agree- 
ment is excellent both in positions of peaks and their 
magnitudes. However, this agreement is lost, if the ten- 
sor force is removed from the cross-shell interaction, as 
shown in the right panel of Fig. [2j For instance, the high- 
est 0c?3/2 peak is shifted in the wrong direction, and the 
main peak of 0d 5 / 2 moves away towards higher energy. 
In the present calculation, as already stated, the ESPEs 
around 40 Ca are consistent with experiments with a rea- 
sonably large 0d 3 / 2 -0G? 5 / 2 gap ~7 MeV. The proton shell 
structure evolves from 40 Ca to 48 Ca, giving rise to the 
agreement with the fragmentation of spectroscopic fac- 
tors. In particular, because only the tensor force can 
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FIG. 2: (Color online) Spectroscopic factors of proton hole 
states measured by 4S Ca(e,e'p) [j| (upper) and its theoret- 
ical calculation (lower left). The cross-shell tensor force is 
removed in lower right panel. The black, blue and red bars 
correspond to lsi/2, 0^3/2 and lds/2 states, respectively. 
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FIG. 3: (Color online) Intuitive illustration of the struc- 
ture of intrinsic state at Z=14. q implies intrinsic quadrupole 
moment (fm 2 ). 



change the 0d 3 /2-0ri5/2 gap to this order of magnitude 
(by ^2 MeV), the agreement shown in Fig. [2] provides 
us with the first evidence from electron scatteringexper- 
iments to the tensor-force-driven shell evolution |2| . 

We now consider shape transitions driven by shell evo- 
lution, by taking as an example the exotic Si (Z=14) 
isotopes with even iV=22~28. In a conventional view, 
Z=14 is a (sub-)magic number with a large 1si/ 2 -0g?5/ 2 
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gap, as shown in Fig. OUa). For a large gap, six protons 
occupy 0^5/2, forming a closed shell. This should end up 
with a spherical shape particularly for a doubly magic 
nucleus, 42 Si (N=28), similar to that observed in 34 Si 
(N=20). The situation changes when many neutrons oc- 
cupy the O/7/9 orbital. Following the mechanism of Ot- 
suka et al. [2J, the monopole interaction of the tensor 
force is strongly repulsive between a proton in 0d 5 i 2 and 
a neutron 0/7/2- Hence, as neutrons occupy O/7/2, the 
spin-orbit splitting decreases and the energy of proton 
0^5/2 l cve l comes closer to lsx/2 level (The ls 1 / 2 energy 
is unaffected by this mechanism). 

Figure [3^b) presents an intuitive illustration on the 
consequence of this change by taking a simple case com- 
prised of lsi/2 and 0g?5/2 orbits. These orbits are mixed, 
like Nilsson orbits, due to the quadrupolc deformation 
of the intrinsic state. Assuming an axially symmetric 
deformation, single-particle states of the same magnetic 
quantum numbers, denoted m, are mixed in the intrin- 
sic states. Figure [3] (b) indicates that this occurs for 
m = ±1/2 between lsi/2 and 0d$/ 2 , with amplitudes 
sinf? and cos#, respectively. A smaller ls 1 / 2 -0d 5 / 2 gap 
results in more mixing, while the phase of mixing ampli- 
tude depends on as to whether the shape is prolate or 
oblate. In the present case, protons occupy the states 
of m = ±5/2, ±3/2, which yield in total a negative in- 
trinsic quadrupolc moment (i.e., an oblate shape). The 
total intrinsic quadrupole moment gains a larger magni- 
tude, if the 1si/2-0g?5/2 mixing gives a negative moment. 
The proton-neutron interaction, apart from its monopole 
part, can be modeled by a quadrupole-quadrupole in- 
teraction, and a similar mixing occurs for neutrons in 
lp 3 / 2 and O/7/2, producing a negative intrinsic moment. 
Thus, by having the mixing leading to a negative in- 
trinsic moment, the total magnitude of the moment 
becomes larger for both protons and neutrons, giving 
rise to stronger binding of the intrinsic state from the 
quadrupole-quadrupole interaction. The actual structure 
is determined by the competition between actual (effec- 
tive) gaps and effects of the quadrupole-quadrupole inter- 
action, implying that the shapes of Si isotopes depends 
on how the gaps are reduced by the tensor force. 

We then examine this picture in the context of the 
shell-model calculations with the SDPF-MU Hamilto- 
nian. Figure 2] exhibits yrast properties of even- A Si and 
S isotopes. Effective charges are (e p ,e n ) = (1.20e, 0.45e) 
fixed already by properties of lighter isotopes. The over- 
all agreement to experiment is remarkable. For instance, 
in the present result, 2+ levels of Si isotopes keep com- 
ing down as N increases consistently with experiment, 
whereas some increase is seen at TV = 28 in other shell- 
model calculations [2(| 27 1. The nice agreement suggests 



that the intuitive picture above holds particularly to- 
wards 7V=28, resulting in strongly deformed shapes with 
low excitation energies consistent also with recent mea- 
surement by GANIL 12J. In fact, if the tensor force 
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FIG. 4: (Color online) (a,b) 2+ and 4+ energy levels and 
(c,d) B(E2;0^2+) values of Si and S isotopes for 7V=22~28. 
Symbols are experimental data [9l-fl5l[. Solid (dashed) lines 
are calculations with (without) the cross-shell tensor force. 



is omitted from the cross-shell interaction, 2+ level of 
42 Si28 goes up. It is of much interest to see the missing 
experimental data in Fig.|4]as well as more precise B(E2) 
values. Figure H] exhibits results for S isotopes in good 
overall agreement. 

The potential energy surface (PES) is useful to under- 
stand shapes contained in theoretical calculations. Fig- 
ure [5] exhibits PES for Si isotopes obtained by the con- 
strained Hartrec-Fock method [Hj for the SDPF-MU 
Hamiltonian. The full Hamiltonian is taken in Fig. [5] 
(a~d), whereas the cross-shell tensor force is removed 
in Fig. [5] (e~h). Shape evolutions are seen very clearly 
in both sequences (a~d) and (e~h), starting with sim- 
ilar patterns in 36 Si. The shape evolves as more neu- 
trons occupy p/-shell, with distinct differences between 
the two sequences. In (b,c), the deformation becomes 
stronger from (a) with triaxial minima, whereas the 
shape becomes more like modestly prolate in (f,g). In 
(d), one finds a strongly oblate shape with a sharp min- 
imum, but the minimum is at the spherical shape in 
(h). This strong oblate deformation produces low 2 + 
level and large B(E2) in Fig. 0] for the "doubly-closed" 
42 Si. Thus, the shape of exotic Si isotopes changes signif- 
icantly within the range of AN ~6. This is considered to 
be a manifestation of Jahn- Teller effect with varying shell 
structure driven by the tensor force. We note that the 
two sequences produce rather similar levels and B(E2) 
in Fig. [4] for lighter isotopes, and the structure of 42 Si 
serves a key role in the study of the tensor- force effect. 

Although more experiments arc needed to clarify 
present issues, there are some hints of the triaxiality 
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FIG. 5: (Color online) Potential energy surfaces of Si isotopes 
from N = 22 to 28 calculated with (left) and without (right) 
the cross-shell tensor force. The energy minima are indicated 
by red circles. 
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level. The present calculations with (without) 
the cross-shell tensor force locates it at 0.62 (1.11) MeV 
above 2+ for 40 Si. The extraordinary low-lying 2^~ level 
caused by the tensor force appears to agree with a recent 
it is proposed that at least either 
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7-ray experiment 
of the observed 7-rays 638(8) and 845(6) keV directly 
feeds the 2+ state. 

With density- functional methods (3fj| - [33j . Si and S are 
discussed with somewhat less overall agreement with ex- 
periment, compared to the present calculation. It will be 
of interest to see single-particle properties given by these 
works in the view of Jahn- Teller effect. For instance, 
how much does the lsi/2-0d 5 / 2 gap change from 40 Ca to 
48 Ca? 

In summary, we have presented more evidences of the 
tensor-force driven shell evolution, proposed by Otsuka 
et al. Q , in low- lying states of K isotopes and in distri- 
bution of spectroscopic factors measured by 48 Ca(e,e'p) 
experiment. Similarly, the levels and B(E2)'s of exotic Si 
and S isotopes are described well by the same Hamilto- 
nian. The nuclear deformation at low excitation energy 
is a Jahn- Teller effect, which should be sensitive to the 
shell structure. The shell evolution driven by the ten- 
sor force plays a crucial role in rapid shape transitions, 
including a robust mechanism for the appearance of the 



stable oblate shape at subshcll closures against prolate 
dominance of the nuclear deformation. In future, simi- 
larities and differences in 78 Ni will be of great interest. 
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